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Abstract. - The iterated random walk is a random process in which a random walker moves 
on a one-dimensional random walk which is itself taking place on a one- dimensional random 
walk, and so on. This process is investigated in the continuum limit using the method of 
moments. When the number of iterations n — > oo, a time-independent asymptotic density is 
obtained. It has a simple symmetric exponential form which is stable against the modification 
of a finite number of iterations. When n is large, the deviation from the stationary density is 
exponentially small in n. The continuum results are compared to Monte Carlo data for the 
discrete iterated random walk. 



Introduction. When a walker moves at random along a straight support, its typical 
displacement at time t grows like / ~ t 1/>2 . If the support of the walk is itself a random walk 
(RW), at time t the walker is typically located at the Zth step of the walk on which it is 
moving. Thus the mean-square displacement of the walker behaves as (X 2 (t)) ~ I ~ t 1 / 2 and 
his actual typical displacement is reduced, growing as t 1 ^ 4 instead of t 1 / 2 for the conventional 
RW. 

This model of a RW on a RW has been studied for a one-dimensional support pQ and also 
in higher dimensions It has been used to discuss the diffusion of the stored length along 

a polymer chain entangled in a network and the Brownian motion of charged particles 

in a turbulent plasma |7j. In this latter case, the particles are constrained to diffuse along the 
random magnetic field lines in the limit B — > oo. 

In this work we study the properties of the one-dimensional iterated random walk (IRW) , 
i.e., the RW on a RW which is itself a RW on a RW, and so on. The first three steps of 
this iteration process are shown in fig. ^ Working in the continuum (Gaussian) limit, we 
first study the moments (X p (t)) of the probability density P^ n '(X,t) to find the walker at 
{X, t) after the nth iteration. In the next section, we examine the limit n — > oo for which 
the asymptotic probability density p(°°\X,t) can be obtained explicitly. We also look at 
the way this asymptotic density is approached when n is large. Then the continuum result 
is compared to Monte Carlo data for the discrete IRW. The universality of the asymptotic 
probability density and possible generalizations are discussed in the final section. 
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Fig. 1 - Construction of the position X of the random walker at time t for the IRW: The walker in (a) 
starts at xo(0) = X(0) =0 and moves at random along the sequence of steps of the RW in (b). Thus 
the coordinate xo(t) is also the time coordinate ti for the support in (b). In the same way, xi gives 
the number of steps performed by the RW in (c) on which the RW in (b) takes place. Finally, 
at tz = X2, the last RW in (d) is located at X, which gives the actual position X(t) of the random 
walker. 



Moments of the iterated probability density. - We assume that our random walker starts 
from X — at t — 0. All other one-dimensional RWs, on which he is moving, extend in time 
from — oo to +oo and have a common origin. In the continuum limit, to the symmetric RW 
is associated the Gaussian density 

p(x,t) = -jL=exp(-?-\ , t > 0, (1) 



IWt V 2t 

and the probabilty density to find the walker at (X, t) after n + 1 iterations is obtained by 
integrating the product of the successive probability densities p(xi, \ti\) with ti = Xi-i, to = t 
for the random walker and x n+ i = X for the last support, over all the intermediate coordinates 
Xi for i = to i = n (see fig. [Q: 

/>+oc /"+OC r + oo 

P (n+1) {X,t) = / dx n p(X, \x n \) / dx n - 1 p(x n ,\x n ^ 1 \) ■ ■ ■ dx p(x 1: \x Q \)p(x ,t) . 



(2) 

As a consequence, the iterated probability density satisfies the recursion relation: 

r+oo 

P^ n+1 \X 1 t)= / dxp(X,\x\)P^(x,t). (3) 



In the following we leave the initial probability density p(°\x, t) unspecified. Since the Gaus- 
sian density p( X 7 \ x\ ) in © is even in X, the same is true of all the iterated densities P( n '(x, t), 
n > 0. Their odd moments vanish while their even moments are given by 



M 2 ( ;»(f)= / dXX 2 vp( n \X,t) = M%>{t), M^{t)= / dX\X\ p P^(X,t), (4) 
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where the moments M.p of \X\ have been introduced. According to ©, these moments may 
be written as 

M p n+l) (t) = dxP( n \x,t) dX \X\ p p(X, \x\) = 2 p l 2 K2 fl , 2 ' / dx\x\ p / 2 P^(x,t), 



where the last relation follows from the form of for the Gaussion density, 



(5) 



dX \X\p )=JU = -^jjf- (2tf/ 2 . (6) 

-00 \f lut 1 [2) 



Equation © leads to the recursion relation 



^(n + l) (t) = 2P /2£i| + l) A1 (n) (t) 



(7) 



so that 

H-'W-^WII^ rn, ~ 2) ■ (8) 

fc=i V2J 

Making use of the identity [8] 

n(,) = fr^±I)= ^i±^ (9) 

M r (3) 2 2 -(i-i/2")r + 1) ' l9j 

the moments of the iterated probability density may be written as: 

Mil? (t) = — -^— 9 M { "\ 9n (t) , M^ ] +1 (t) = . (10) 

Asymptotic probability density and how it is approached. - Since the walker and all the 
underlying RWs have a common origin in space and time, 

P M [X,0) = S(X). (11) 

When t > and n — > 00, the even moments in (llUfl considerably simplify, 

M^(t>0) = M, (12) 

and become time-independent since M { °\t) = 1 for any initial density function. The Fourier 
transform of the asymptotic probability density is related to its moments through a Taylor 
expansion: 

dXe ikX P {oc ' ) {x,t> 0) = J2 ^-M^°°\t > 0) 



p=0 
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The inverse Fourier transform, 

i r+ao 2 f+oo —ikX 

P^(X,t>0) = — dke- ikx V (co) {k,t>0) = - dk — -, (14) 

is easily evaluated using the method of residues. Finally, the asymtotic density takes the 
simple form 

P<~>(*,«)={ e S W III , (15) 

which is time-independent as soon as t > and universal, i.e., independent of the form of 
the initial density p(°\ This asymptotic density is invariant under the Gaussian transforma- 
tion ©. 

Let us now examine how this stationary density is approached when n is large. For this, 
the initial density P^ must be specified. Coming back to the random-walk problem, we use 
the Gaussian density for which 

M°u*) = £ T[fr {2t)p/2 " (16) 

according to The even moments in I|1U|) take the following form: 

M ln )(t) (2p)! 2*/»T(£ + $) /2 „ 



Using the duplication formula [S], 

r(2*) _ o22 - 1 r(z + |) 
r(z) r(i) 

this expression can be rewritten as: 



2p 1 j " 4p r (i + £) 



(18) 



(19) 



Thus Mj™ (0) = 5 p .o in agreement with (|11|) and a length scales as t Y l 2 + after n iterations. 
When n 3> 1, a first-order expansion in powers of 2~" gives 

(2£) p/2 " = l+ln(2t)^+0(4-") 

r(i + |r) = i + *(i)| r + 0(4-«) = i- 7 A + o(4-») 

^(t) = ^{l + [ln(2i)+ 7 ]^}+0(4-»), (20) 

where ^(z) = T'(z) /r(z) is the digamma function and 7 = .577 215 664 ... is Euler's constant. 
Thus, the leading correction to the Fourier transform of the asymptotic density reads 

«><«>(*,*> o) = in(2 : )+ ^f> f4y = - in(2 : )+7 j fc2 ^ 2 (2D 

v y 2™ ^ V 4 / 2™ (fc 2 +4) 2 v ; 
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Fig. 2 - Scaling function for the iterated probability distributions of the discrete RW compared to 
the stationary density (solid line) obtained in the continuum limit. The histograms of the probability 
distributions have been obtained through Monte Carlo simulations for t = 10 2 to 10 5 and t+1, working 
with 10 7 samples. The statistical fluctuations are much smaller than the symbols as indicated by the 
symmetry of the data. 



and the inverse Fourier transform leads to: 

6PM(X,t> 0) = QX \ - 1) e-^l . (22) 

The deviation from the asymptotic density is exponentially small in n, a behaviour which is 
independent of the initial density function according to <|1U[I . It has a weak (logarithmic) time 
dependence for the Gaussian density. 

Scaling behaviour and discrete IRW. - The Gaussian density in Q can be written under 
the more general scaling form 

_ ,-l/z - ' X 



p ( Xl t)=t-y*f{— z ) , (23) 

where the dynamical exponent z, which is equal to 2 for the Gaussian density, governs the 
transformation of the time, t' — t/b z , under a change of the length scale, x' = x/b. The 
scaling behaviour of the iterated density can be obtained by induction. Assuming that 



(24) 



l/z 

which is true with zq = z and n — for the initial density, and making use of J2J), we obtain 
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so that p( n+1 ) scales like P^ in (|24p. The dynamical exponent evolves according to z n — 
zz n -i — z n+1 and the scaling function satisfies the recursion relation 

Thus, with z > 1, the asymptotic density p(°°}(X,t) is time-independent and reduces to its 
scaling function f(°°'(X). 

We now make use of this scaling behaviour to analyse Monte Carlo data for the discrete 
IRW. The continuum result is expected to give a good description of these data for not too 
small values of X and t. The histograms for the discrete probability distributions have been 
obtained as follows: First a RW with t steps is generated, leading the walker to Xo (see fig-DJ. 
Then a second walk with \xq\ steps, corresponding to the first support, is constructed. When 
it ends at x\ the value of N^(xi), giving the number of walks ending at x% at the first 
iteration, is updated and the same process is repeated for the following iterations. N s = 10 7 
samples have been generated in this way for t = 10 2 , 10 3 , 10 4 , 10 5 . Since all the x n have the 
parity of t, the simulations were repeated for t + 1 in order to double the number of points 
along the X axis. The properly normalized histograms for a comparison with the density 
functions are obtained by dividing by 2N S . 

The scaling functions for the discrete iterated walk are compared to the invariant den- 
sity ()15[) in fig. • A good data collapse is obtained for each value of n. As expected, there 
is a rapid convergence to the asymptotic density. The convergence is from above when \u\ is 
greater than a value close to 1, in agreement with l|22l) . 

Final remarks. - We have seen that the Gaussian iteration process leads to an invari- 
ant probability density l(15jl which is independent of the initial density P' '. Actually, the 
universality of the invariant density is more extended since P^ may be itself considered as 
resulting from a finite number of iterations involving arbitrary densities pk(x,t) (k — 0, to). 
The initial density being now po(xo,t), we have: 

r+oc r+oc 

P (0) (ir,£) = / dx m ^ip m (x, |x m _i|) / dx m -2Pm-i(x m -i, |aJm-a|) • • • 



-foe 

dx a p 1 (x 1 , \x \)p a (x ,t) . (27) 

-oo 

After an infinite number of Gaussian iterations this density will evolve, as before, to the 
universal time-independent invariant density. The arbitrary densities may be, for example, 
shifted Gaussian densities 



/ s 1 

Pk(x,t) = — ==exp 



(x -v k tf 



2t 



(28) 



corresponding to RWs with arbitrary drifts v^. 

When, in the iteration process, the Gaussian density is replaced by a density which is 
even in x, scales like in (|23|) and has finite moments, the nonvanishing even moments of the 
iterated density are easily obtained in the same way as above. They involve a product similar 
to that appearing in (JHJ and are time-independent when n — > oo. Unfortunately, we were 
unable to find another instance where these moments could be simplified as in 110(1 . allowing 
for an explicit calculation of the time-independent asymptotic density. 
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